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We consider diffusion processes with the help of Markov random walk models. Especially the 
process of diffusion of a relativistic particle in a relativistic equilibrium system is considered. We 
interpret one of the results as causal Zeno effect for its similarity to quantum Zeno effect. Another 
problem we considered is about collision number. Basing on our numerical results, we propose that 
in the considered situation the probability density distribution among different collision numbers is 
a lognormal distribution. 



I. INTRODUCTION 

Diffusion process is very common in physics, chemistry, 
biology and many other fields. It can be described by a 
diffusion equation with or without a finite maximum ve- 
locity P,i,ll,[3. Another approach toward diffusion pro- 
cess is numerical simulation [1, 0, Yet, it can also be 
considered with the help of random walk models (RWM) . 
Many literatures take this way, such as [1, [^, • 

In this letter, we consider the diffusion of a par- 
ticle in an equilibrium system in the framework of 
Markov random walk model. We begin our study 
with the onc-particlc velocity distribution of a non- 
relativistic/relativistic equilibrium system. 

As it is well-known, the onc-particlc velocity distribu- 
tion of a ncar-frcc non-rclativistic equilibrium gas is the 
Maxwell's probability density function (PDF): 

/M(v;m,/3) = (f^)''/^e-^-^/^ (1) 

ZTT 

The first relativistic generalization of Maxwell's function 
was proposed by F. Jiittner from a consideration of max- 
imum entropy principle [111, [H, : 

/jp(p;™,/3) = ^e-^^, (2) 

where E = ^ + = ?7i7(v), p ^ mv7(v), and 
7(v) = l/Vl — v^. Zj is the normalization factor. Writ- 
ten in velocity PDF, Jiittner's distribution is 

/7.(v;m,/3) = ^7(v)^+''e-'3"-^(-). (3) 
Zj 

A lot of work has been done basing on Jiittner's dis- 
tribution since its proposal. However, in the 80's of 
the last century, doubt about Jiittner's function was ex- 
pressed and "modified Jiittner's function" was proposed 
[isl . \\A . [Til l. But recent numerical simulations favored 



Jiittner distribution Q as the correct one-particle veloc- 
ity distribution of a relativistic equilibrium system 0, 0] . 

So, in this letter Maxwell's function ([1]) and Jiittner's 
function ([3]) are used. Numerical analysis is done. By 
comparing the results from the two distributions, we in- 
terpret one of the results as causal Zeno effect for its 
similarity with the Quantum Zeno Effect [l^, . 

Another problem we consider is that how many times 
a particle will collide with other particles during the pro- 
cess of diffusion from a initial point to a final point in a 
given time interval. As one can conceive, there will be 
a probability distribution among different collision num- 
bers. One of our results is that if the final point is not 
specified, the probabilities to collide whatever N times 
will be the same. For the case of that the final position is 
also fixed, we make numerical analysis and propose that 
the probability distribution is the lognormal distribution 



II. RANDOM WALK MODELS 

In this section, we consider diffusions as random walk 
processes. At first we give some general descriptions, 
then we apply the model to the non-rclativistic Markov 
diffusions, finally we consider the relativistic Markov dif- 
fusions. Our discussions will be confined to the one di- 
mensional case. 



A. General descriptions 

When a particle is diffusing in an environment, one 
wants to know the probability for it to reach some posi- 
tion from a given position within a given time interval. 
The formal kinematics formula is 

i(t) = xo + / dsv{s), (4) 
Jo 
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where v[^s) is its velocity. 

For a free particle, its velocity is a constant. While if 
it interacts with other matters, its velocity will change. 
Here, we assume that the particle interacts with other 
particles only by point collision. Thus the above formula 



2 



transforms into 



N 



(5) 



i=l 



where ti,i = 1,2, N — 1 arc the tunes of collisions and 
In is the given ending time. Wc write Ti = ti — for 
simplicity. Although the actual collision number is not 
N but A'^ — 1, one can roughly say it is or it is indicated 
by iV. 

{{vi\; {Ti\) defines a path in the space. If the particle 
has a probability density /({«,}; {ti]) to follow this path. 
Then the transition probability density function of order 
N can be written as 



(6) 



PN{t,x\Q,XQ) — J dvi...dvN J dTi...dTN 

fi{v^h{nmx~m)■ 



If the collision number is not concerned about, the total 
transition PDF is 



p{t, x\0, Xq) = ^ pN{t, x\0, Xo). 



(7) 



In the following, we set ti be equal (= t = t/N). And 
we consider Markov process in which case f{{vi};{Ti}) 

can be written as YliLi /("fi! '''i)- 
Thus the N-order transition PDF is 



N 



6{x-m (8) 



PN{t,x\0,xa) = 
Using 

S{x) = — / dke''''' 
together with (O, one can get 

-1 /"OO r ^ r- 

{t,x\Q,xo) = — y dfce*'=(^-^«)^n j dvJ{vf,T) 



(9) 



Pn 



If the velocity probability density function of Vi is fi (w) . 
then the PDF of Vsn is 

fsNiV) = hiv) * h{v) * ... * fNiv), (13) 

here * means convolution. 

This leads to the probability density function of Vn' 

hiNiV) = fsN{NV)N. (14) 

Finally the N-order transition PDF is 

X — Xq 1 

PN{t,x\0,XQ) = fMN{ )-• (15) 

We will follow this way when we make numerical anal- 
ysis. 



B. Non-relativistic Markov diffusions 

Firstly we apply this model to the case of non- 
relativistic Markov diffusions. The equilibrium velocity 
PDF of a non-relativistic system is the Maxwell distribu- 
tion. 

Substituting ^ into pH]) . and assuming 



/3 = 



1 



ksT 2mD ' 



(16) 



one gets 



PN{t, x\0,xo) 



1 



-g 4DNt 



y/AnDNr 

1 (x-jjq)^ 



(17) 



2ti 



dfce*'=(^-^«V(-'tT) 



N 



(10) 



where ^p{u) is the characteristic function of f{v), i.e. its 
Fourier transformation. 

One can also get an expression for pN{t,x\Q,xo) by 
invoking the mean velocity. Define 



N 



Vsn = ^Vi, 



(11) 



which is the solution to the ordinary diffusion equation 
Jl S with a proper initial condition. 

One can easily see that this expression of transition 
J'DF is in conflict with the special relativity, for when 
|a; — a;o| > ct, there is still a small but non- vanishing 
probability for the particle to be found. 

Another property of non-relativistic Markov diffusions 
is that for different values of N, the pjv(t, xjO, Xo)'s are 
the same, i.e. no matter how much the collision number 
is, even when it is infinite, pN{t,x\Q,xa) is the same as 
above. So, the particle can always diffuse. As we will 
see, this is not the case in relativistic Markov diffusions. 



Relativistic Markov diffusions 



then the mean velocity 

_ 1 ^ 1 ^ 1 

Vn^-^Y. ^''^^ = ^ E = M^SN. (12) 



1. Formula development 

For relativistic Markov difi^usions, one should use 
Jiittner's distribution ^ instead of Maxwell's distribu- 



3 



tion. Then the N-order transition PDF will be 



PNit,x\0,xo) 
1 

1 

2^ 



dvfj{v]T)e 



-ikrv 



N 



N 



(18) 



This is very difficult, or perhaps, even impossible to 
carry out the integral analytically. For a compensation, 
we consider its large-N asymptotic behavior by invoking 
its N-order mean velocity V n- 

The N-order mean velocity during the time interval 
from to t is 



y N — — 



1 ^ 



1 ^ 

N ^ 



Vi. 



(19) 



i=l 



Since M{vi) = 0, yar(vi) = < oo. Central Limit 
Theorem (CLT) [1^ asserts that the distribution of 



"^^^iVi/aVN converges to Gauss's normal distribution 

(20) 

(21) 



with parameters (0,1): 

Nix;0,l) = 

Then, the probability density about Vn will be 

1 ~w 



^27rcr2 /N 
And the N-order transition PDF will be 

1 _(x-3if 



PN{t,x\0,xa) 



v27rfT(7^ 



(22) 



Here is the variance of Jiittner' distribution (|Aip : 
it is a function of x(= "m-P), and t = t/N. 

At a given reciprocal temperature (3, is fixed. When 
iV ^ oo, T ^ 0. Then 



PAr(t, a;|0, xo) 6{x~ xq). 



(23) 



This result means that the particle will not diffuse at 
all in this situation. By comparing the results in non- 
relativistic and relativistic situations, one can see that 
they are contrary to each other. In the non-relativistic 
situation, we have seen that pN{t, x\0, xq) is independent 
of N, and always is a Gaussian distribution. But here, in 
the relativistic case, we have the result (l23l). 



"freeze" the evolution of the system by measuring it fre- 
quently enough in its (known) initial state [l^ . 

In the random walk model, N ^ oo means that the 
particle is continuously collided by other particles. In 
other words, it is continuously observed. So the result 
(|23p means that the particle will never depart from its 
initial position if it is continuously observed. However, as 
we have seen, in the non-relativistic situation, the particle 
can always diffuse. 

This difference between the two situations arises from 
the different velocity distribution functions. Further- 
more, as one can see, when x becomes large, Jiittner's 
distribution ([3]) tends to Maxwell's distribution, save for 
its finite support. Thus the boundedness of velocity is 
indispensable for the convergence to a delta function. 

One can also perceive that at a given temperature, 
there will be a probability distribution among different 
values of N, and the behavior of pN{t, x\Q,xo) with the 
most possible value of N will dominate in the diffusion 
process of the particle. 



3. Xn = 1 

As for the probability distribution among different val- 
ues of N, we first consider the following expression: 



A 



TV 



PNit, x\0, XQ)dx. 



(24) 



This expression quantifies the probability of the particle 
colliding iV — 1 times with other particles in a given time 
interval t — Q from the initial position xq to non-specified 
final positions. 

Substituting (fTU|) into the above expression, the out- 
come is 



Xn — 



1 1'°° 

dx— / dke'^^''-''«'^Lpj(-kTf . (25) 



27r , 



As 



one gets 



6(k) = — / dxe*'=("-"«\ 



2. Causal Zeno effect 

We venture to interpret the result (|23p as causal Zeno 
effect for its similarity with the quantum Zeno effect. 
Quantum Zeno effect is a name coined by George Su- 
darshan and Baidyanath Mishra in 1977 in their analy- 
sis of the situation in which an unstable particle, if ob- 
served continuously, will never decay . One can nearly 



which means that the probabilities for the particle to 
collide different TV — 1 times with other particles are the 
same. And it is right even for cases that t; (i=l,2,...,N) 
are different and the process is non-Markovian. 

It is surprising as we have thought that there would 
be a non-trivial probability distribution. However, in the 
case with a specified final position, one may anticipate 
the probabilities will be different. We will check this by 
numerical analysis. 
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FIG. 1: p]v(i, a;|0, 0) against x with x=0-l s^nd t=l. This is 
the ultra-relativistic case. 
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FIG. 2: pM{t, x\0,0) against x with x=10- We have used 
Juttner's function. This is the non-relativistic case. However 
the resuhs are stiU different from (|17|l which is independent 
of N. Here large N sharpens the peak because the support of 
fjiv) is finite. 
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FIG. 3; The probabihty distribution among different values of 
N before normalization. The three curves are of x=0.1, 0.08, 
0.05, respectively. One can see that these curves are similar 
to the curves of the lognormal distributions(See Fig|6}. 



III. NUMERICAL ANALYSIS 

A. Pjv(t, a;|0, 2:0) against x 

The formulae we use for numerical calculation have 
been presented in Section [TTl In the following, we use 
Juttner's function and we set xo=0- 

For ppf{t,x\0,0) agaist x, there are three parameters: 
X for the velocity distribution; t for time of the diffusing 
process; N for the kind of the diffusing pathos. 

We give our results in Fig.([T]) and Fig.©. In Fig.([T]), 
the value of x is 0.1. The S peak in the middle belongs to 
the curve of pn with N=2. This curve has another two 
peaks at ±1. The two second highest peaks which are 
symmetrically located belong to the curve of N=3. This 
curve has low sharp peaks at ±1 too. For N=4, there 
are 3 peaks symmetrically located besides the peaks at 
±1. This case is not shown in the picture. For large N, 
the pN{t-, x\0,0) arc Gaussian-like, which is predicted by 
CLT. 

In Fig. ([2]), X — 10- This is the case of low- velocity. 
The peaks are more sharpened. The effective velocity of 
the diffusing particle is small and the particle is more 
confined in its neighborhood. And with N increasing 
the peak sharpens. This is different with the result pT]) 
which comes from Maxwell's function for the PDF of ve- 
locities. Jiittner's function with large x applies to the 
low velocity situations as well as Maxwell's function, but 
since the support of Jiittner's function is finite, the causal 
structure is maintained. 



B. Pisi{t, x\0, xq) against N 

For pN{t, x\0, 0) against N, the parameters are x, t and 

X. 

In the ultra-relativistic situation and in the region with 
X far away from zero, the distribution among different 
values of N will be characterized by sharp peaks. Because 
in the region with x far away from zero, pN{t, x\0,0) is 
characterized by sharp peaks (Sec Fig.([T])). If the pa- 
rameter X is fallen into a peak of some pN{t^ x\0,0), the 
probability distribution to this value of N will be very 
large, and those to others will be very small. 

As for the case with x near zero, contributions to 
the total transition PDF (the sum of pAr(t, ccjO, 0)'s with 
different values of N) come mostly from large N. How- 
ever, for large N, p]y{t, x\0,0) is gaussian-like. So we 
just consider the large x cases in which pAr(t, x|0, 0)'s are 
Gaussian-like. 

We give the results in Fig. ([3]). From Fig.Q, one can 
see that the curves in this picture are much similar to the 
curves of the lognormal distributions in Fig. ([6|) . Lognor- 
mal distribution has been widely used in chemistry, bi- 
ology, ecolo gy, social sciences and economics, and many 
other fields [191 . We propose that the probability distri- 
bution among different values of N is a lognormal distri- 
bution. 
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FIG. 4: Jiittner's function. When x is big, the temperature is 
small and fj{v) reduces to a Gaussian-like distribution. For 
small X the peaks appear at large velocities. 
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FIG. 6: Lognormal distribution 
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FIG. 5: The variance and the standard deviation of Jiittner's 
distribution, i.e. and a, respectively, against the parameter 



APPENDIX A: JUTTNER'S DISTRIBUTION 

Jiittner's distribution in 1-dimension space (see Fig.Q 
and Fig.®) is 

/j.(«;m,/5) = ^^^7(«)'e->^"^''\ (Al) 

where x = /^'ti, Kiix) is a modified Besscl function of 
the second kind. 

The mean is and the variance is 

~ / d?;7(i;)3e-x'>'(") 



IV. CONCLUSIONS 

In this letter, we have considered random walk models 
of diffusion processes. One of our results is that there 
will be a causal Zeno effect in causal diffusions. As for 
the problem of collision number, we have found that the 
probabilities for the diffusing particle to collide whatever 
A'' — 1 time with other particles in a given time interval 
will be the same if the initial position is specified but the 
final position is not. For the case of that the final position 
is also fixed, we have made numerical analysis and pro- 
posed that the probability distribution among different 
values of N is a lognormal distribution. 
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APPENDIX B: LOGNORMAL DISTRIBUTION 

Let (7^,o-2 (x) be Gauss's normal distribution, 

1 _ (x-AOf. 



V27rCT2 



then 



f(^\ _ / h9^l,<7^-il'09x), if X > 



if X < 0. 



(Bl) 



(B2) 



is called the Lognormal distribution with parameters 
li,a'^{-oo < yu < oo, < < oo) [H, [ll, [131 (see 
Fig.®). 
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